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IN-SITU  INVESTIGATION  OF  TEMPERATURE  AND  BIAS  DEPENDENT 
EFFECTS  ON  THE  OXIDE  GROWTH  OF  Si  AND  Ge  IN  AN  ECR  PLASMA 
Y.Z.  Hu,  Y.Q.  Wang,  M.  Li,  J.  Joseph  and  E.A.  Irene 
Department  of  Chemistry,  CB#  3290 
University  of  North  Carolina 
Chapel  Hill,  NC  27599-3290 

ABSTRACT 

The  electron  cyclotron  resonance,  ECR,  plasma  oxidation  of  Si  and  Ge  was  investigated  using 
in-situ  spectroscopic  ellipsometry  at  substrate  temperatures  of  80  to  400  °C  and  at  bias  voltages 
of  -30  to  +60  V.  A  study  of  the  oxide  growth  kinetics  by  ECR  plasma  oxidation  results  in  three 
distinct  regions  of  g-owth  with  the  first  two  being  linear  and  the  last  parabolic.  At  the  earliest 
linear  stage  the  rate  of  oxide  growth  is  the  fastest,  and  corresponds  to  around  3  nm  film 
thickness  which  is  not  dependent  on  bias.  Following  this,  the  second  linear  region  displays  an 
oxide  growth  rate  proportional  to  the  bias  with  typical  growth  rates  of  0.10,  0.32  and  0.60 
nm/min  for  0,  +30  and  +60  V,  respectively,  at  300°C.  The  third  region  displays  parabolic 
kinetics  and  corresponds  to  the  Cabrera-Mott,  C-M,  theory  for  the  oxidation  by  charged  species 
in  the  limit  of  a  low  electric  field.  Activation  energies  of  0. 19  and  0.28  eV  are  obtained  using 
the  C-M  model  for  the  ECR  plasma  oxidation  of  Si  and  Ge,  respectively.  Increasing  the  negative 

bias  causes  oxide  sputter  etching  and  surface  damage,  with  no  significant  contribution  to  ?or 
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oxidation.  A  probable  mechanism  for  the  ECR  plasma  oxidation  process  is  also  discussed.  O 
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I.  INTRODUCTION 


Low  temperature  processing  of  semiconductors  is  of  interest  since  rates  for  thermally 
activated  defect  production  and  redistribution  of  impurities  are  greatly  reduced.  Plasma 
oxidation,  utilizing  an  oxygen  plasma  is  a  promising  low  temperature  technique  used  to  grow 
dielectric  films  on  semiconductor  surfaces.  Microwave  electron  cyclotron  resonance, ECR, 
plasmas  have  much  higher  fractional  ionization(up  to  10%),  longer  ion  and  energetic  neutral 
mean  free  paths,  and  lower  maximum  ion  energies  than  conventional  rf  parallel-plate 
discharges1.  Recently  an  ECR  oxygen  plasma  has  been  used  to  grow  silicon  dioxide  from  single 
crystal  silicon2'5,  exhibiting  the  physical,  chemical  and  electrical  properties  similar  to  those  of 
high  temperature  thermally  grown  oxides.  However,  reports  on  the  ECR  plasma  oxidation  of  Ge 
appear  to  be  lacking.  Although  the  ECR  plasma  oxidation  of  semiconductor  materials  are  of 
great  interest,  the  kinetics  and  mechanisms  involved  are  complicated  and  thus  are  not  yet  fully 
understood. 

The  work  presented  in  this  paper  focuses  on  a  kinetic  study  of  the  various  parameters 
influencing  the  oxidation  process  of  Si  and  Ge,  using  both  in-situ  static  spectroscopic 
ellipsometry,  SE,  and  dynamic  real  time  single  wavelength  ellipsometry,  SWE.  The  dependence 
of  oxide  growth  on  time,  substrate  temperature  and  bias  voltage  are  presented  along  with  a 
growth  model  based  on  that  of  Cabrera  and  Mott5,6. 

II.  EXPERIMENTAL  TECHNIQUE 

The  ECR  plasma  oxidation  equipment  used  in  this  study  was  described  previously5.  A  2.45 
GHz  microwave,  as  generated  by  a  magnetron,  is  guided  through  a  quartz  window  to  the  ECR 


plasma  source  chamber  filled  with  oxygen  gas  (99.997%)  that  was  taken  directly  from  a  cylinder 
without  further  purification.  The  vacuum  chamber  can  be  evacuated  to  5xl0'8  Torr  using  a 
turbomolecular  pump.  The  substrate  to  be  oxidized  is  placed  on  a  stainless  steel  holder  that  can 
be  heated  using  a  halogen  lamp.  The  temperature  was  monitored  and  controlled  using  a 
thermocouple  attached  to  the  front  side  of  the  sample  holder  and  adjacent  to  the  sample.  A 
calibration  of  the  holder  temperature  was  made  in  comparison  with  substrate  temperature.  The 
oxygen  plasma  irradiates  the  substrate  at  near  normal  incidence. 

The  plasma  oxidation  conditions  were  as  follows:  total  pressure  of  lxlO'3  Torr,  oxygen  flow¬ 
rate  of  20  seem,  microwave  power  of  300  W,  substrate  temperature  from  25°C  to  400°C,  and 
substrate  bias  voltages  from  -60  V  to  4-60  V.  The  samples  used  were  commercially  available 
lightly  doped  P-type  <  100 >  Ge  wafers  and  P-type  <  100 >  Si  wafers. 

The  thicknesses  of  the  oxide  layers  was  determined  by  means  of  real  time  single  wavelength 
ellipsometry.  The  various  thicknesses  reported  in  this  work  were  obtained  by  a  trajectory 
method5,7.  The  single  wavelengths  chosen  for  Si  and  Ge  are  340  and  335  nm,  since  the 
measurements  are  insensitive  to  temperature  and  highly  sensitive  to  oxide  growth  at  these 
wavelengths.  The  SE  measurements  were  taken  at  41  photon  energies  between  2.5  and  4.5  eV. 
For  the  measurement  of  each  sample  the  offsets  of  analyzer  and  polarizer  were  redetermined. 
This  was  necessary  to  maintain  the  highest  accuracy,  because  a  small  misalignment  of  the  sample 
stage  occurred  with  sample  changes.  Each  spectroscopic  measurement  takes  about  10  min  and 
during  the  measurement  the  plasma  was  stopped.  For  the  oxides  grown  above  room  temperature 
the  spectra  were  obtained  at  room  temperature  after  cooling  the  sample.  This  was  done  because 
of  the  unavailability  of  reliable  optical  data  for  all  the  components  present  for  the  higher 
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temperatures. 


HI.  RESULTS  AND  DISCUSSION 
A.  Optical  Modeling 

Spectroscopic  ellipsometry  has  been  used  extensively  to  study  thin  films  grown  on 
semiconductor  materials8"12.  The  ellipsometric  measurement  yields  two  angles:  A  and  Sk,  from 
which  the  complex  reflection  ratio  pcxp  is  obtained:  p„p=tanS^eiA.  From  a  literature  database  for 
the  known  constituents  of  the  film  and  substrate,  peti  is  calculated  and  compared  with  p„p  and 
as  figure  of  merit  for  comparison,  an  unbiased  estimator  <5  is  calculated  from  the  relationship513: 


6  = 


N 

X/  IPexp'Pcail 


1/2 


U) 


here  N  is  the  number  of  wavelengths  sampled,  and  P  the  number  of  unknown  parameters.  A 
minimizing  procedure  gives  the  best  fit  parameters  which  are  film  thicknesses  and  volume 
percents  for  the  constituents  at  the  90%  confidence  level. 

Figure  1  shows  a  typical  variation  of  the  pseudo-dielectric  function,  <e> ,  as  obtained  from 
SE,  in  terms  of  both  s{  and  e2  versus  photon  energy,  for  the  plasma  oxidation  of  Si  and  Ge. 
Four  sensible  models  shown  in  Fig. 2  were  evaluated,  and  model  D  made  up  of  pure  oxide  as 
the  topmost  layer  and  an  interface  layer  composed  of  oxide  and  amorphous  Si  or  Ge,  yields  the 
best  fit  in  terms  of  the  lowest  5.  The  calculation  results  from  model  D  for  both  Si  and  Ge  are 
also  shown  in  Fig.l.  It  is  interesting  to  observe  that  the  best  fit  model  for  both  Si  and  Ge  ECR 
oxidation  is  the  same. 
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B.  Time  Dependence  of  the  Oxide  Growth 

In  order  to  investigate  the  ECR  plasma  oxidation  kinetics,  the  plasma  oxidation  of  Si  and 
Ge  was  performed  for  various  times  at  different  bias  voltages  and  temperatures  under  the  same 
plasma  conditions.  Figure  3  shows  the  typical  results  for  the  oxide  thicknesses  as  a  function  of 
time  at  a  microwave  power  of  300  W,  substrate  temperature  of  300  °C,  and  three  different 
substrate  bias  voltages  of  0,  +30  and  +60  V.  The  results  of  Fig. 3a  and  3c  show  the  evolution 
of  the  thickness  as  a  function  of  the  oxidation  time  ,  obtained  by  the  trajectory  method  and  a  one 
film  model  from  the  real  time  single  wavelength  ellipsometric  measurements  during  the  first  20 
and  30  min  for  Si  and  Ge,  respectively.  The  thickness  results  of  Fig. 3b  and  3d  were  obtained 
by  SE  using  a  two  film  model.  The  total  thickness  is  given  by  the  solid  line  and  the  interface 
thickness  by  the  filled  symbols.  The  experimental  results  in  Fig. 3  show  that  there  are  three 
regimes  for  the  relationship  between  time  and  thickness.  The  earliest  two  regimes  are  nearly 
linear  with  the  first  being  both  the  fastest  and  virtually  bias  independent.  Following  these  linear 
regimes,  a  parabolic  shape  in  the  growth  data  is  observed. 

In  order  to  discuss  the  time  dependence  of  the  oxide  growth,  several  oxidation  models, 
including  the  linear-parabolic  model  for  Si  thermal  oxidation  developed  by  Deal  and  Grove14  (D- 
G),  a  time  power  law  oxidation  model  by  Wolters  and  Zegers-van  Duynhoven15  (W-Z)  and  the 
parabolic  model  by  Cabrera-Mott616  (C-M),  have  been  used  J-51718.  However,  we  were  not 
successful  with  any  of  above  models  for  the  entire  experimental  thickness-time  range  of  our 
plasma  oxidation  data. 

Figure  4  shows  a  schematic  of  the  DC  potential  profiles  in  the  growing  oxide  as  a  function 
of  the  thickness3.  Initially,  the  surface  potential  V,  will  be  equal  to  the  bias  potential  V„,  i.e., 
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for  bare  wafers  or  for  oxides  that  are  too  thin  to  sustain  an  electric  field.  After  the  oxide  has 
reached  a  thickness  Lox  the  overall  oxide  potential  drop  Vox=  Vb-V,  will  increase  with  increasing 
film  thickness.  Therefore,  the  field  in  the  oxide  will  be:  E0I=V0I/L0I.  Previous  experimental 
results  also  indicated  that  negatively  charged  oxygen  atoms  formed  at  the  plasma-oxide  interface 
are  the  primary  oxidant  under  anodic  oxidation  conditions2,3,5,12,17.  The  field  Eox  may  lead  to  rapid 
transport  of  ionic  species  through  the  oxide. 

Since  the  oxidizing  species  react  at  the  plasma-wafer  interface  without  transport  in  the  first 
linear  regime  where  Eox=0,  the  oxidation  is  bias  independent.  In  the  second  linear  regime  the 
electric  field  develops  in  the  thin  oxide  layer,  and  the  oxidizing  ionic  species  transports  rapidly 
by  the  field  due  to  low  transport  resistance.  Thus  the  film  growth  remains  reaction  limited,  but 
it  now  is  electric  field  dependent.  These  two  linear  regimes  in  the  early  ECR  plasma  oxidation 
of  Si  and  Ge  are  very  similar  to  those  formed  in  the  early  stage  of  Si  thermal  oxidation19,20. 
From  the  theory  by  Murali  and  Murarka20  (M-M)  at  the  earliest  stage,  the  oxidizing  species 
diffuse  through  a  very  thin  oxide  layer  into  the  substrate  due  to  the  low  diffusion  resistance  of 
the  thin  oxide  layer.  The  oxidation  takes  place  over  a  reaction  volume  zone  rather  than  only  at 
the  interface,  and  the  rate  of  oxide  growth  is  fast.  When  the  oxide  thickens,  the  oxidizing  species 
is  consumed  at  the  interface  due  to  an  increased  diffusion  resistance  of  a  thicker  oxide  layer. 
According  to  M-M  model  there  is  a  critical  thickness,  Lox’,  below  which  the  oxide  growth  is 
given  by  a  high  growth  rate.  For  the  thicknesses  greater  that  than  Lox*  the  growth  is  given  by 
a  second  slower  rate.  The  critical  thickness  is  not  bias  dependent,  but  it  is  dependent  on  the 
substrate  temperature.  From  Fig.  3a  and  3c  and  our  previous  paper5  typical  values  of  Lox’  are 
about  1.5  nm  (for  Si)  at  room  temperature  and  about  3  nm  (for  Si  and  Ge)  at  300  °C.  Following 
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these  linear  regimes,  parabolic  growth  is  observed  which  agrees  with  the  C-M  model5  6 16  which 
includes  electric  field  effects  resulting  from  the  applied  biases. 

C.  Bias  Dependence  of  the  Oxide  Growth 

The  reactive  species  generated  in  oxygen  plasmas  are  positive  ions,  negative  ions,  electrons 
and  neutrals.  Positive  or  negative  bias  applied  on  the  substrate  causes  different  effects,  since  the 
effected  reactive  species  is  different.  For  a  constant  positive  applied  bias  voltage,  the  dominant 
reactive  species  are  likely  to  be  O  and  O. 

According  to  the  C-M  model  for  oxidation  by  charged  species  in  the  limit  of  a  low  electric 
field,  the  growth  rate  is  given  by  the  relationship: 


Lox  ~  C.exp 


'-'a 

KT 


Vnr 

+  C2 
KT  2 


(2) 


where  Vox  is  the  potential  drop  across  the  oxide  film,  E,  is  the  thermal  activation  energy 
associated  with  the  diffusivity  of  oxidizing  species  in  the  applied  electric  field,  t  the  time,  K  the 
Boltzman  constant,  T  the  absolute  temperature,  and  C]  and  C2  are  constants.  In  order  to  compare 
our  experimental  results  to  C-M  model,  the  dependence  of  the  parabolic  coefficient,  B,  with  the 
bias,  shown  in  Fig. 5b  and  5d,  is  examined  and  is  given  as: 


B  = 


d(L 


dt 


2, 

ox  '  =  k  V  +  C 

Kb  Vb+Lb 


(3) 


where  Vb  is  bias  voltage,  Cb  and  k,,  are  constants.  Figure  6  shows  the  parabolic  coefficient 
dependence  to  be  approximately  linear  with  bias  and  using  the  values  of  ^=0.25  nm2/min*V 
for  Ge  and  1^=0.023  nm3/min*V  for  Si.  It  is  clear  from  eq.(2)  and  eq.(3)  that  the  relationship 
between  Vox  and  the  bias,  Vb,  is  approximately  linear.  The  increase  of  the  growth  rate  with  the 
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positive  bias  (on  the  Ge  and  Si)  indicates  that  the  negative  oxygen  ions  are  the  principal  oxidant 
species.  These  ions  may  come  from  the  plasma  itself  or  be  formed  at  the  surface  via  electron 
attachment  to  the  neutral  02  and  O.  The  attachment  is  a  favored  process  and  causes 
decomposition  of  Oj.  Evidence  for  this  has  been  observed  for  thermal  oxidation21  and  electron 
stimulated  oxidation  of  Si22,23.  At  negative  bias  positive  oxygen  ions  and  atoms  are  the  likely 
dominant  reactive  species.  An  increase  of  the  negative  bias  results  in  oxide  sputter  etching  and 
surface  damage. 

D.  Temperature  Dependence  of  the  Oxide  Growth 

The  temperature  dependence  of  the  oxide  film  growth  is  determined  from  the  data  in  Fig. 5a 
and  5c  at  oxidation  temperatures  from  room  temperature  to  400°C.  The  C-M  model  predicts  that 
the  parabolic  coefficient  follows  an  Arrhenius  expression,  and  this  is  confirmed  for  our  results 
in  Fig. 7a.  The  linear  fit  of  these  data  yields  an  apparent  activation  energy,  E,.  The  values  are 
0.28  eV  and  0.19  eV  for  the  plasma  oxidation  of  Ge  and  Si,  respectively.  It  is  interesting  to 
notice  that  the  oxide  growth  rate  on  Si  is  higher  than  that  for  Ge  at  low  temperatures  for  ECR 
plasma  oxidation.  It  is  well  known  that  the  oxide  growth  rate  of  Si  is  much  smaller  than  Ge  in 
thermal  oxidation.  For  a  comparison  between  plasma  and  thermal  oxidation,  we  measured  the 
Ge  oxide  growth  rate  at  the  oxidation  temperatures  from  400  to  550  °C.  Figure  7b  shows  the 
parabolic  coefficient  dependence  on  temperature  for  the  thermal  oxidation  of  Si  and  Ge,  with 
the  value  for  Si  at  about  one  order  of  magnitude  smaller  than  that  of  Ge.  This  clearly  indicates 
a  very  different  oxidation  mechanism  for  the  ECR  and  thermal  processes.  The  FCR  plasma 
environment  contains  atomic  charged  oxidant  species,  and  as  indicated  above  in  the  bias  effect 
measurements  the  negative  oxidant  species,  likely  O'  are  causative.  In  the  fields  existing  in  the 
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ECR  system  under  bias  and  no  bias  conditions  these  charged  species  are  affected  by  the  electric 
field.  Thus,  the  transport  of  the  smaller  O"  species,  as  compared  with  molecular  0:  for  thermal 
oxidation,  in  both  a  concentration  and  field  gradient  is  probably  the  reason  for  the  large 
difference  in  the  activation  energies. 

W.  CONCLUSIONS 

A  comparison  of  the  ECR  plasma  oxidation  between  Ge  and  Si  shows  that  the  ECR  plasma 
oxidation  is  very  similar  to  thermal  oxidation  in  some  points  and  greatly  different  for  other 
points. 

The  similarities  are: 

1.  The  oxide  structure  is  two  layer:  a  pure  oxide  on  top  and  an  interface  layer. 

2.  The  kinetics  of  oxide  growth  shows  three  oxidation  regimes:  two  linear  and  one  parabolic. 

3.  The  oxide  properties  are  similar. 

The  differences  are: 

1.  The  dominant  oxidizing  species  is  negative  oxygen  ions  for  plasma,  but  molecules  for  thermal 
oxidation. 

2.  The  oxidation  rate  is  bias  dependent. 

3.  The  activation  energies  in  the  plasma  oxidation  of  Si  and  Ge  are  much  smaller  than  those  in 
thermal  oxidation. 

4.  In  low  temperature  ECR  plasma  oxidation,  the  oxide  growth  rate  of  Si  is  higher  than  that  for 


Ge  which  is  reversed  for  thermal  oxidation. 
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FIGURE  CAPTIONS 


Figure  1.  Pseudo-dielectric  function  of  Si  (a)  and  Ge  (b).  Experimental  data  after  ECR  plasma 
oxidation  are  presented  by  symbols.  The  solid  lines  are  before  oxidation.  The  dashed  lines 
correspond  to  an  two  layer  model  (model  D  in  Fig. 2)  calculation. 

Figure  2.  Four  optical  models  used  to  evaluate  the  SE  data  of  Fig.  1 .  The  top  is  for  Si  oxidation 
and  the  bottom  is  for  Ge  oxidation. 

Figure  3.  Thicknesses  of  the  grown  oxide  layer  vs  ECR  plasma  oxidation  time  for  three  biases, 
(a)  and  (b)  are  for  Si  oxidation  at  room  temperature;  (c)  and  (d)  are  for  Ge  oxidation  at  300  °C; 
(a)  and  (c)  Thicknesses  are  calculated  at  single  wavelength  using  a  one  layer  model;  (b)  and  (d) 
Thicknesses  are  calculated  from  SE  data  using  model  D  in  Fig.2.  The  solid  lines  are  for  total 
thicknesses.  The  interface  layer  thicknesses  are  represented  by  filled  symbols. 

Figure  4.  Schematic  diagram  of  potential  profiles  in  the  growing  oxide  at  constant  positive  b:as 
voltage:  (a)  Initial  growth;  (b)  thicker  oxide. 

Figure  5.  Square  of  total  oxide  thickness  vs  the  oxidation  time  from  SE  data  for  various 
temperatures  and  biases:  (a)  Si  at  +30  V;  (b)  Si  at  room  temperature;  (c)  Ge  at  +20  V;  (d)  Ge 
at  300  °C. 

Figure  6.  Parabolic  coefficients  deduced  from  linear  fit  of  the  data  in  Fig. 5(b)  and  (d)  vs  applied 
bias:  (a)  Si;  (b)  Ge. 

Figure  7.  Arrhenius  plots  of  the  parabolic  coefficients:  (a)  plasma  oxidation  of  Si  and  Ge,  data 
deduced  from  Fig. 5  (a)  and  (c);  (b)  thermal  oxidation  of  Si  and  Ge.  The  data  for  Si  are  from 
Ref.  14. 
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